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Abstract 

Spin-spin interaction causes the mixing between ground state wave functions of 
baryons containing three quarks of different flavors. We examine the effect of this 
i-£h 1 mixing on the baryon masses in the framework of the modified bag model. 
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Invented more than 30 years ago the MIT bag model [Tf2f3] still serves as an 
useful method which yields reasonable predictions for a variety of hadronic 
properties (at least for the ground states). Various aspects of the Bag model 
are discussed in a number of review papers [41I51I6] . Originally the model was de- 
signed for the ultrarelativistic case of the light quarks and was rather success- 
ful in describing the low-lying hadron spectrum. Although the first straight- 
forward application of the model to calculate the spectrum of the hadrons 
containing heavy quarks was of very limited success because the disagreement 
between calculated and observed data was evident, later on the bag model was 
adjusted to incorporate the heavy quarks. It was recognized that the reconcil- 
iation between the bag model and the heavy quark physics could be achieved 
by taking into account the so-called c.m.m. (center of mass motion) correction. 
For the hadrons containing only one heavy quark there exist an approximate 
solution of the c.m.m. problem - one can simply associate the center of mass 
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with the heavy quark and fix it 08]. For the baryons containing two heavy 
quarks one can proceed in a similar way assuming the heavy quarks to form 
the doubly- heavy diquark and then put this object at the center of the bag 
[9]. The advantage of such approach is the simple and clear physical picture, 
but the price for this simplicity is three more ore less related bag models: 
one for the hadrons consisting of light quarks (in this case the role of c.m.m. 
correction is partly played by the so-called zero-point energy), one for hadrons 
containing single heavy quark, and one for the baryons containing two heavy 
quarks. Moreover, there remains the case of the baryons consisting of three 
heavy quarks which needs special treatment. Another rather popular method 
to deal with c.m.m. problem is to employ a wave-packet ansatz [10]. Both 
approaches give only approximate solutions to the problem. However, the sec- 
ond seems to be more universal and could be preferable in the case one tries 
to obtain the unified description of the light and heavy hadrons. The c.m.m. 
corrections is not all the story, and in order to have plausible unified descrip- 
tion of light and heavy hadrons (mesons and baryons) in the framework of the 
bag model some other QCD inspired improvements such as running coupling 
constant and running quark masses are necessary [TT] . 

When the number of quark flavors increases we are confronted with additional 
problem which needs some clarification. For the spin-1/2 baryons containing 
three quarks of different flavors there exist two states with the same spin and 
parity. We can construct the set of orthogonal wave functions by assuming 
the first two quarks to be in the relative spin-0 or spin-1 state, respectively. 
In general case the physical states would be the linear combinations of these 
mathematical states 



The mixing mechanism depends on the model and approximation used. In the 
MIT bag it would be the hyperfine color- magnetic interaction. The same is 
true in the ordinary constituent quark model as well as in almost all variants 
of the potential model. Some estimate of the state mixing is also possible in 
the heavy quark effective theory [T2] . 

Because of the ambiguity in how the quarks are to be ordered in the math- 
ematical wave functions \(qiq 2 ) s q 3 ) the expansion in Eqs. ([T]) is not unique. 
If one does not want to bother about the quark ordering, one can simply di- 
agonalize the interaction energy matrix |13| . Since such a procedure has not 
become a common practice, the very natural question arises, how some au- 
thors have managed to avoid this state mixing problem. The answer was given 
a long time ago in the Ref . [2] . The authors of that paper have shown that for 



\B)=C 1 \(q 1 q 2 )°q 3 ) + C 2 \(q 1 q 2 ) 1 q 3 ), 



By = C 2 \(q 1 q 2 ) q 3 )-C 1 \(q 1 q 2 ) 1 q 3 ). 
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the interaction energies with quark mass dependence ~ l/{m a mb) an optimal 
quark ordering scheme could be found. The prescription is to pick the closest 
in mass quarks as the first two in the wave function | ((71(72) 93)- Then the mix- 
ing of the states with different S values is small, and the effect of this mixing 
on baryon masses is negligible. Strictly speaking, such mixing exists even in 
the light baryon sector between the wave functions from which the physical 
states S° and A are constructed. Because of the approximate isotopic sym- 
metry this mixing is small, and if we are interested only in the calculation of 
baryon masses, we can safely ignore this effect. The explicit calculations with 
the isospin-symmetry violating terms taken into account show that the mixing 
is indeed small |15j, as expected. In the sector of charmed baryons there are 
also rather strong indications that the mixing between S c and "E' c baryons is 
small, with negligible shifts in the masses of these hadrons again [IB] . 

And what could be said about the wave function mixing in the framework 
of the bag model? Of course, we expect that all reliable models of hadron 
structure yield similar results. However, we cannot apply the results of Ref. [Ti] 
directly because the dependence of the interaction energy on quark masses 
in the bag model is somewhat more complicated. For example, the values 
of light quark masses in the bag model could be set to zero, while in the 
nonrelativistic models these values approach one-third of the nucleon mass. 
A simple way to make the things clear is to perform direct calculations in 
the bag model taking the mixing interaction into account. This means that 
in the calculations of baryon energy the off-diagonal matrix elements of the 
color-magnetic interaction should be included. 

In this paper we are going to examine the mixing of the ground state baryon 
wave functions in the framework of the modified MIT bag model. In the next 
section we give a short description of the model we are dealing with. The 
concluding section contains the results of our investigation accompanied by 
the discussion and some additional remarks on the validity of approaches with 
and without mixing. 



2 The model 



The ground state energy of the hadron defined in the static spherical cavity 
approximation is given by 

Asjr 

E = —BR 3 + Y,n l s l + AE, (2) 

where B is the bag constant, R is the bag radius, Si is eigenenergy of the 2th 
quark in the cavity, and AE stands for the interaction energy. AE consists of 
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Table 1 

Parameters which specify the color-magnetic interaction energy of baryons consist- 
ing of three distinct quarks. 
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color-electric and color-magnetic parts as described, for example, in Ref. [TT] 
in detail. For our purpose the most important is the contribution of the color- 
magnetic interaction, which in the case of the baryons containing three distinct 
quarks can be written as 

AE m = a c (R) aijMij (m h rrij , R) . (3) 

j>i 



Here a c (R) is the running strong coupling constant. The function My (m^ rrij, R) 
depends on quark masses and hadron bag radius and it can be calculated 
explicitly. Parameters dy- specify the spin dependence of the interaction en- 
ergy between quarks qi and qj. They are proportional to the matrix elements 
((.lilz) 81 qzli.o'i - & j)\(qiq2) S2 qs) i where <7j are appropriate spin generators. These 
coefficients can be calculated straightforwardly using algebraic technique, as 
described in the Ref. [T7|, and the transformation of the basis [IB] 



l(<M 2 ) Jl2 g 3 > J = E (-l)^ 3+Jl2+Jl3 y/(2J 12 + 1) (2J 13 + 1) (4) 

Jis 

x < } |(gig 3 ) Jl3 g 2 ) J , 

{ 33 J Jl3 J 

where necessary. The results are presented in the Table [TJ where for simplicity 
the abbreviations \ J\2) = |(<?i<?2) Jl2 <73) are used. 

The relation between the calculated bag-model energy E and the hadron mass 
M is given by 

E = J d i s<l> 2 P {s)VM 2 + s 2 , (5) 
where $p(s) is a Gauss profile 
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The effective momentum P specifies the momentum distribution and is defined 



as 



P 2 = lJ2^Pl (7) 



Here pi are the momenta of the quarks. The c.m.m. parameter 7 is to be 
determined in the fitting procedure. For the baryons containing 6-quarks the 
relation (jSj) gives practically the same results as the familiar Einstein relation 

m 

M 2 = E 2 — P 2 . (8) 



In the presence of 6-quarks we prefer to use this simple relation instead of 
rather cumbersome procedure based on Eq. (15"]) . For the running coupling 
constant a c (R) and running quark mass rfif(R) we use the following expres- 
sions: 

rfif(R) =m f + a c (R) -5f , (10) 



where Rq is the scale parameter analogous to QCD constant A. Parameter A 
helps us to avoid divergences when R — ► R Q . For each quark flavor we have 
two free parameters m/ and 5/ to be adjusted. 

Let as summarize our zoo of model parameters. These are the bag constant 
B, the c.m.m. parameter 7 which determines the strength of the c.m.m. cor- 
rections, two parameters (A and Rq) from the running coupling constant 
parametrization, and finally six parameters im s , 5 S , m c , 5 C , m,},, 5^) necessary 
to define the running quark mass functions rrTf(R). The light (up and down) 
quarks are taken to be massless. To fix the parameters B, 7, A, and Rq the 
experimentally observed masses of the light hadrons (N, A, tt, and the average 
mass of the uj—p system) were chosen. To fix the mass function parameters rhf, 
Sf for each quark flavor we have employed the masses of corresponding lightest 
vector mesons {(f), J/ip, T) and the mass values of the lightest baryons Af con- 
taining the quark qj of the corresponding flavor. We employ the same fitting 
procedure as in our previous work [11] , and the values of the parameters to be 
used as the input in the bag model calculations are: B = 7.597 • 10~ 4 GeV 4 , 
R = 2.543 GeV~\ A = 1.070, 7 = 1.958, m s = 0.161 GeV, S s = 0.156 GeV, 
m c = 1.458 GeV, 5 C = 0.112 GeV, m b = 4.793 GeV, 5 b = 0.061 GeV. The 
parameters B, R , A, 7, m s , 6 S are the same as in Ref. [llj. The numerical 
values of the remaining four parameters (m c , S c , mb, Sf,) differ slightly from the 
corresponding values presented in [TT] because in the present work we have 
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Table 2 

Dependence of the calculated energies before the matrix diagonalization (in GeV) 
on the arangement of quarks for the S c , E/ c (columns 2-4) and S b , r! b (columns 5-7) 
baryons. The last two rows contain the squared expansion coefficients Cf, Cf of the 
wave functions obtained after matrix diagonalization. 



(<7l<72)<?3 
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(sc)u 
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(ub)s 


(sb)u 


E(l) 
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6.013 


6.010 


E(0) 
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6.059 


6.062 


cl 


0.9950 


0.3139 


0.1912 
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0.2656 


0.2347 


cl 


0.0050 


0.6861 


0.8088 


0.0003 


0.7344 


0.7653 



used new more accurate values of A c (2.286 GeV) [20] and A b (5.620 GeV) 
[2T] masses. 



3 Results and discussion 

Let us proceed to the discussion of our main point of concern - the wave 
function mixing of heavy baryons in the bag model calculations. The ground 
state baryons we are interested in are E c , E' c , E*; E b , E' b , E b ; E bc , E bc , S£ c ; 
and Q bc , Q! bc , f^ c . The mixing is possible only between the spin-1/2 states (H c 
and E' c , for example). In order to calculate the masses of all these baryons we 
use the bag model parameters listed at the end of the preceding section. For 
the spin-3/2 states denoted as | . . .)* the calculation procedure is exactly the 
same as adopted in the paper [TT]. We minimize the energy Eb* of each such 
baryon as a function of the bag radius R and then apply Eq. ([5]) (for E*) or 
Eq. (jSJ) (for El, S£ c , and f^ c ) to determine the corresponding baryon masses. 
For the spin-1/2 states | . . .) and | . . .)' the procedure differs only in the choice 
of the energy function to be minimized. In this case we use the trace of the 
energy matrix Eb + Eb> which remains invariant under state mixing. Then 
we calculate the diagonal and off-diagonal matrix elements of the interaction 
energy, diagonalize the energy matrix, and use Eq. <$J5§ or (JSJ) again to deter- 
mine the masses of the physical baryons. To gain some insight how the things 
look like we present some intermediate results of calculations in the Tables [21 
[3j In the first two rows of these tables we give the c.m.m. uncorrected energy 
values E(l) and E(0) corresponding to the mathematical wave functions in 
which the first two quarks in the spin coupling scheme {qiq-i) qz are in the 
spin-1 and spin-0 states. The last two rows contain the squared wave function 
expansion coefficients obtained after matrix diagonalization. The symbols b, 
c, s denote the bottom, charmed, and strange quarks, respectively, and for the 
sake of simplicity the symbol u is used for both light (up or down) quarks. 
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Table 3 

Dependence of the calculated energies before the matrix diagonalization (in GeV) 
on the arangement of quarks for the Eft c , E' bc (columns 2-4) and Ofe c , Q' bc (columns 
5-7) baryons. The last two rows contain the squared expansion coefficients Cf, Cf 
of the wave functions obtained after matrix diagonalization. 



(<7i<72)<73 


(uc)b 


(ub)c 


(cb)u 


(sc)b 


(sb)c 


{cb)s 


E(l) 


7.078 


7.050 


7.041 
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7.224 
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E(0) 


7.035 


7.062 


7.072 


7.212 


7.235 


7.241 


cl 


0.9833 


0.3693 


0.1474 


0.9872 


0.3538 


0.1591 


cl 


0.0167 


0.6307 


0.8526 


0.0128 


0.6462 


0.8409 



The inspection of results in Tables |2] and [3] shows a striking dependence of the 
calculated energies on the quark ordering. As one can see, the wave function 
with two first quarks in the relative spin-0 state has the lowermost energy 
only when the heaviest quark (e.g., 6-quark) is picked up as the third in the 
corresponding spin coupling scheme {q\q2) S Qz- This is the only case when the 
traditional prescription 

|...> = I(<M 2 ) <?3>, |...> , = |(9ife) 1 fe>. (11) 

could be maintained, because the energy matrix diagonalization leads to neg- 
ligible changes of the initial energy values. The direct calculations show that 
even in the most problematic case of the E bc — E' bc system the difference be- 
tween energy values before and after diagonalization does not exceed 1 MeV 
and is obviously much smaller than the systematic uncertainties of the model. 
So, if one is interested only in the baryon mass spectra, one can adopt the 
prescription fill I) , construct the optimal basis by arranging the quarks in in- 
creasing order of their masses, and not to bother about the diagonalization of 
the energy matrix anymore. At this point a remark is necessary. One must be 
very cautious when dealing with other baryon parameters (such as magnetic 
moments, for example). As it was shown in Ref. [Hj, the wave function mixing 
may change the values of the calculated magnetic moments substantially even 
when the optimal basis is used. Although this problem is beyond the scope 
of the present paper, it is worth attention, and we are going to return to this 
question in the future. 

Before going to the concluding remarks we want to compare the masses of 
baryons calculated in our work with the results obtained in other models and 
experimental data where available. We have chosen for the sake of compari- 
son the baryon mass estimates in nonrelativistic [22j and relativistic [231124^ 
potential models obtained in the quark-diquark approximation, the estimates 
obtained in the quark-diquark approximation of the bag model [9], the cal- 
culations in the simplified variational approach [25], and predictions provided 
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Table 4 

Masses of S c , E' c , E* and S&, H£ baryons (in GeV). The row denoted as Bag 
contains the results obtained in our work. The row Exp contains averaged over the 
isodoublet experimental energy values. 
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m 


2.474 


2.578 


2.655 


5.808 


5.946 
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[26|27] 


2.468 


2.582 


2.651 


5.810 


5.955 


5.984 


Exp 


2.469 


2.577 


2.646 









Table 5 

Masses of E bc , Hj , Hj and f2{, c , ^{, c > ^tc baryons (in GeV). The row denoted as Bag 
contains the results obtained in our work. 



Particle 


-be 


1=1 be 


^bc 


^bc 




"fee 


Bag 


6.846 


6.891 


6.919 


6.999 


7.036 


7.063 


[22] 


6.82 


6.85 


6.90 


6.93 


6.97 


7.00 


m 


6.933 


6.963 


6.980 


7.088 


7.116 


7.130 


[26] 


7.029 


7.053 


7.083 


7.126 


7.148 


7.165 


M 


6.838 


7.028 


6.989 


6.941 


7.116 


7.077 



using various sum rules based partially on the heavy quark symmetry con- 
siderations [26.27J. The experimental values are taken from the Particle Data 
Tables [28]. The data for the baryons of Eq type are presented in the Table H] 
and for the Eq x q 2 , flQ 1 Q 2 type baryons in the Table [5j 

From Table H] it is seen that for the baryons containing one heavy quark all 
approaches give rather similar qualitative picture. Inspection of the Eq — E*q 
hyperfme mass splitting indicates that in our version of the bag model the 
interaction energies for these baryons could be slightly underestimated. Com- 
parison with experiment also shows that all approaches give reasonable results. 
One could even insist that owing to the approximate nature of the models the 
agreement with experiment (though not excellent) is surprisingly good. Such 
success gives us some confidence that we are on the right path in understand- 
ing the properties of heavy baryons. 

For the baryons with two heavy quarks the situation is somewhat different. 
As seen from Table El all but one approaches give similar qualitative pictures 
of the baryon spectra again. A striking exception is the results obtained in 
the paper [9] (the reversed order of the E' bc , El c and Q' bc , f^ c states). The bag 
model results for the ground state baryon masses calculated in our work are 
laid out somewhat above the estimates [22] obtained in the nonrelativistic 
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potential model based on the quark-diquark approximation. Relativistic ap- 
proach [23] gives similar mass spectrum as ours but shifted approximately 70 
MeV upwards. The predictions based on the sum rules are higher than our 
estimates by approximately 170 MeV and 110 MeV for the S;, c and f2& c families 
respectively. The difference between the baryon spectrum obtained in the pa- 
per |9J and the others is of qualitative character. It could look strange, but it is 
the direct consequence of the attempts to incorporate the mixing effects for the 
ground state baryons in the quark-diquark approximation to the bag model. 
We already know that in the ordinary approach the wave function mixing can 
play an important role in the calculations of the baryon mass spectra. However, 
in general, we cannot draw a direct link between the quark ordering in the 
spin coupling scheme and the corresponding diquark structure. Nevertheless, 
some correspondence between the two pictures is expected. For example, in 
the quark-diquark approximation to the potential model the physical E' bc and 
fl' bc states are those with scalar cb diquark [24J , as could be expected from the 
analogy with the ordinary approach (see the 4th and 7th columns of Table [3]). 
In the usual approach the interaction of the system consisting of two quarks 
with the third one is provided by the interaction of its individual constituents. 
The mixing of the wave functions is possible only when the interaction be- 
tween the first and the third quarks is of different strength as compared with 
the interaction between the second and the third (as the quark becomes heav- 
ier its hyperfine interaction with other quarks decreases). When the mixing is 
present the correct mass splitting is achieved only after the diagonalization of 
the energy matrix. On the other hand, in the quark-diquark approximation 
some information about the initial structure of the diquark is lost, and, as 
a rule, the mixing of the ground state functions is absent (23]. Maybe some 
remnant of the mixing interaction of the ground states could exist, but prac- 
tically it seems to be not necessary. Since the baryon masses predicted in the 
paper [9J differ radically from the predictions obtained in the bag model with 
the state mixing effects taken into account (this work) and from the results 
obtained in the quark-diquark approximations to the potential model, it seems 
that the mixing effects in the work [9] have been heavily overestimated. Of 
course, in the calculation of energies of the excited baryons one is confronted 
with the mixing of various states, and in consistent calculations [22,24] these 
mixing effects are taken into account. 

As regards the results obtained in our work, first of all, we conclude that, as 
expected, the bag model shares many features of ordinary quark model. The 
main aim of this paper was to examine the heavy baryon ground state wave 
function mixing due to the color-magnetic interaction in the framework of 
the modified bag model. We have found that the main features of the mixing 
interaction in the bag model are the same as in the ordinary nonrelativistic 
quark model. So, we can conclude that fully relativistic treatment of the light 
quarks in the bag has only minor influence on the state mixing properties. For 
the baryons consisting of three quarks of different flavor we cannot in gen- 
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eral ignore the wave function mixing induced by the hyperfine color-magnetic 
interaction. It can even cause sizable changes of the calculated hadronic prop- 
erties. On the other hand, the widely accepted optimal basis can be built up 
by simply choosing the heaviest quark as the third one in the corresponding 
spin coupling scheme. The matrix of the interaction energy in this basis is 
approximately diagonal, and therefore the mixing effects in the baryon mass 
(energy) calculations can be neglected. If for any reason other than optimal 
basis is used, even in the baryon mass calculations the mixing effects must be 
taken into account. 
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